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Here we make use of the Hellmann-Feynman theorem, with the aim to calculate macroscopic 
quantum velocities instead of forces, and we show how it is possible to derive the expression of the 
work per unit time per unit volume (power density) done by a static external uniform electromagnetic 
field interacting with a quantal body. 



INTRODUCTION 

The problem of the dielectric and magnetic response of 
matter, in condensed matter physics, have been inten- 
sively studied in the last decades (and references 
therein) . 

Here we show that making use of a reformulation of the 
Hellmann-Feynman theory, it is possible to derive several 
useful physical properties of a crystalline material in a 
direct and simple way. As a matter of example we shall 
show how to derive the expression of the group velocity of 
an electronic system interacting with an external uniform 
magnetostatic field, and how to obtain the expression of 
the electronic power density. 

Let us start our discussion from a fundamental level 
by briefly reviewing the definition of stationary state or 
equivalently steady state of a quantal body (QB) , as bet- 
ter stated in [Hill Ei. 

Let us consider the time-dependent Schrodinger equa- 
tion 



dt 2m 



V{r,t)^ 



(1) 



where ^ is the state of the system, T the kinetic energy 
operator and V the potential energy. 

A considerable simplification of equation ^ is ob- 
tained if the potential energy is time independent, i.e 
V{r,t) — V{r). In this case it is possible to express 
the general solution of eq.(Q as a sum of products of 
functions of r and t separately. Without entering into 
details, the time dependence of the wave function can be 
expressed as 



*(r, t) = u{r)e ^~ 



(2) 



being u(r ) the eigen function of the following differential 
equation |l8l | 



Eg. (13 formally defines a stationary state (see also 0, 
Mm) for a system where the Hamiltonian assumes the 
form H = T + V. A more general definition of a steady 
state is 



(4) 



where H is the Hamiltonian of a quantal bod y, a nd H 
is either self-adjoint or hermitian (see [l^ lislll^). In 

the following we shall frequently make use of the concept 
of steady states. 



HELLMANN-FEYNMAN THEOREM 

In 1939 R.P. FeynmanE3| developed formulas to cal- 
culate forces in a quantum system, that nowadays take 
the name of Hellmann-Feynman theorem|l5l|. 

The theorem states that for a given configuration of 
nuclei, it is possible to calculate directly the force re- 
quired to hold them; where it is understood that the 
nuclei are to be held fixed in position as point charges. 
Also, given any arbitrary number of parameters A that 
specifies nuclear positions the expression of the force (in 
a generalized sense jl5j) is 



(5) 



where H is the (self-adjoint) Hamiltonian of the sys- 
tem, and ^ are differentiable wave functions. 

Let us here re-state the theorem in order to avoid its 
repetition in the text (see also [isj). 

The wave function must be normalized at every A, i.e. 



(i) / = 1 



(6) 



2m 



V^^- -f V{v)u{v) = Eu{v) 



(3) 



Steady state condition requires 



and E is the corresponding eigen value. 



{ii)H^ = 



(7) 



2 



where E, H and have to be thought depending on A. where dv is a volume element of space, and iJ^ is obtained 
The requirement that _ff is a self-adjoint operator means as 



{Hi) J dv^*HVx^ = J dv\/x-^H^* 
and it is also valid the following relation 



in 



The expectation value, is 



E = / dv^*H^ 



(8) 



(9) 



(10) 



that implies 



VxE = j du«'*VAi?«'+ j dvVx^*H'^+ J du^'*i?VA* 

(11) 

From (i),(ii),(iii), and (iv) we obtain 

VxE = J dv^*VxH^ = J dv<i'*^[inVx,H]^ (12) 

We shall show in the next section, how to make use of 
the theorem (I12II from another perspective. 

GROUP VELOCITY 

In this section it is shown how to calculate the group 
velocity of a quantal body. 

We assume that the differentiable state 4' of the sys- 
tem, varies parametrically with respect to a real vectorial 
field f3. 

We represent ^E* by complex differentiable functions, 
normalized in a suitable region D of space. 

Also, we assume that the Hamiltonian H of the quantal 
body is independent of (3 and self-adjoint. 

Let us define a generalized group velocity vector as: 



Vi3 = n-'VpE{(3) 



We consider here 



where /(/?) is a real vectorial field function of f3. 
E is defined as follows 



(13) 



(14) 



E{l3) "^plHl^p >=< up\Hp\uf} >= I u*pHi3Ufjdv 

(15) 



g-»/(/3)-r^gj/(/3).: 



(16) 



Hp represents the Hamiltonian of the system 30] in the 
Hilbert space spanned by up, and it is parametrically 
dependent on f3. 

Inserting ea. (|15|l in ea. p3|l we obtain, 

nVp ^Vp< ^p\H\^p > = V/3 < uplHplup > (17) 
We may consider Vp to be the averaged vector 



Vp = -< up\VpHp\up > 



(18) 



To prove that, under steady-state conditions, both def- 
initions ca. H18|l and ea H13|) . of group velocity become ex- 
actly equivalent, it is sufficient to make use of ea. (|12|l . 
where we consider here A being substituted by a Carte- 
sian component of p. In the next section, as a matter of 
example, we show implications of ea. H18|) in the theory 
of infinite periodic systems. 



INFINITE PERIODIC SYSTEM 

For an application of ea. (|18|l . to an infinite periodic 
system, it is sufficient to make use of the k — q represen- 
tation, introduced by Zak 3] . 

Here k is the crystal momentum ( 19j) and varies in 
the Brillouin zone of the reciprocal space while the quasi- 
coordinate q varies in the unit cell. 

The k — q representation of Schrodinger's equation, in 
the case of no external fields present in the system, is: 



y(q)]Ck(q)-eCk(q) (19) 



where 



Ho 



2m' 



-ihVc 



(20) 



as also stated in eq.(14) of reference Q, setting explicitly 
to zero the external fields. 

Performing the following phase transformation 



Ck(q) 



wk(q) 



(21) 



the equation for Mk(q) becomes 



[— (-inVq + hk)^ + l/(q)]uk(q) = euk(q) (22) 
Zm 



3 



where 



Ho,u = T^{-in\/^ + hkf + V{q) (23) 

If we now consider a static external electromagnetic 
field {Eo,B^) interacting with the QB, Schrodinger's 
equation in the k — q representation becomes 



2m c 

+ V{q) + eEo ■ iVk]uk(q) 

= eMk(q) 



(24) 



differently from the result shown by Zak, because of 
our choice of the vector potential. 

In fact, in the r representation, we choose the vector 
potential as follows: 



A(r) ^ B" xr 



(25) 



differently from the usual one used for the symmetric 
gauge 26] . 

It is clear that we can apply straightforwardly ea. ((TH|l . 
once we think of the region D as being the super-cell, 
and we substitute f{(3) = /9 = k, and r with q. 

We can also evaluate ea.(l^. at low energies, expand- 
ing brackets (...)^ of eq.||2U in powers of ^. 

Taking only those terms proportional to - , we obtain 



being Gm any r ecip rocal lattice vector and any 
lattice vector[ill9ll2a|. 

Accordingly to boundary conditions on functions 
Uk(q) are 



and 

Uk+G„(q) =e-^^'"'Juk(q) 

for every Gm and Rm- 

A particular solution then reads 



(30) 



Ck(q) 



Uk(q) 



(31) 



Making use of ea. (|18|l and ea. l|24|) . the expression of 
the group velocity operator now becomes 



n 



— X zVk 

rnc 



(32) 



POWER 



Here we show another possible application of ea. H32|l . 
By ea. (l32ll . we can calculate the group velocity of the 
system (in the case of = 0) as follows 



1 P V 

[— (-iWq + nkf + — .{B°x iVk) 

2m c 
+ V{q) + eEo ■ iVk]uk(q) 

= eMk(q) 



(26) 



where v = ■;^[~*^^q + and the Hamiltonian of 
the problem is 



~ 1 cv 
i/k = 7^(-«?iVq + nk)^ + — -{B"x iVk) 
2m c 

+ y(q) + eEo ■ iVk = Ho,), + oj 
The operator a) is instead defined as 



(27) 



d) = eEo ■ iVk 



(B° X iVk) 



(28) 



Boundary conditions on the Bloch-like functions are 
the following 



CkH 



Ck(q + R„ 



Gn 



q) = C^k(q) 

and 

»kI^-Ck(q) 



(29) 



Vnk = < "nk|Vki?k|Unk > 

=< u„k|T^[iVk,i?k]|u«k > 
defining the density current vector as[l9| 



jnk = eVni, 



(33) 



(34) 



We can then define the work, per unit time and per 
unit volume, made by the external electric field Eo on 
the QB as 



(35) 



where fikp is the volume of k points such that their 
associated energies are smaller then the Fermi energy. A 
generalization of eq. H35|) to the case ^ will be given 
in the next section. 
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POWER IN A UNIFORM STATIC 
ELECTROMAGNETIC FIELD 

By ea. (|S^ . noting that the group velocity is not expHc- 
itly dependent on a uniform electrostatic field Eg even 
being eventually present in the system, and bearing in 
mind ea. H34() and ea. (|35|l . we can express the work per 
unit time and per unit volume done by the electromag- 
netic field on the body as 



,E° 



(36) 



Sd 



being S = ;^(Eo x B°) the Poynting vector|33| asso- 
ciated to the electromagnetic field, and 



d = ^ < u„,kNVk|u„,k > 



ed 



dk 
(2^ 



(37) 



(38) 



where e is the electronic charge, and P is the macro- 
scopic polarization^. Ea. l|T7|l and ea. (|SS|l can be eas- 
ily related to the dipole moment per unit volume of 
the systemla, that is expressed as a sum of Berry's 
connections 13 • A review of the physical _and mathemat- 
ical meaning of eq. can be found in 



CONCLUSIONS 

Within our model, concepts of group velocity and 
power density are directly related. Expression H3()|) may 
be useful for calculations of the measured power, ad- 
sorbed or released |U 0| by a quantal body in a 
transformation connecting two of its steady states; it 
may also be useful for calculations of Joule's effect, as 
the latter is directly related to the concept of power, and 
a quantification of the heat acquired or released by a QB, 
along a phase transformation between two of their steady 
states. Ea. H35ll . can be helpful for a more accurate the- 
oretical description of materials that show an hysteresis 
loop'2l'|,i.e during a phase transformation that changes 
their electrostatic potential in analogy with|7i|. As a mat- 
ter of example, we may think to apply our results to the 
study of dielectric properties of materials appertaining to 
the phenomenological polar class |2ljV 
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